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Abstract. - We address the cross effects between mechanical strains and magnetic fields on 
the plastic response of magnetoelastic amorphous solids. It is well known that plasticity in non- 
magnetic amorphous solids under external strain 7 is dominated by the co-dimension 1 saddle-node 
bifurcation in which an eigenvalue of the Hessian matrix vanishes at 7p like y/^p — 7. This square- 
root singularity determines much of the statistical physics of elasto-plasticity, and in particular 
that of the stress-strain curves under athermal-quasistatic conditions. In this Letter we discuss 
the much richer physics that can be expected in magnetic amorphous solids. Firstly, magnetic 
amorphous solids exhibit co-dimension 2 plastic instabilities, when an external strain and an 
external magnetic field are applied simultaneously. Secondly, the phase diagrams promise a rich 
array of new effects that have been barely studied; this opens up a novel and extremely rich 
research program for magnetoplastic materials. 
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Introduction: The well known magnetostriction effect 
where switching on a magnetic field changes the volume 
of a sample of magnetoelastic solid and the Villari effect 
in which mechanical strain changes the magnetization are 
just two examples of the rich variety of cross effects that 
exist in magnetoelastic materials when both an external 
mechanical strains and a magnetic field arc employed. Sur- 
prisingly, the fundamental physics of plasticity has been 
studied much more extensively in the context of pure me- 
chanical strains in non-magnetic amorphous solids, with 
many experiments and numerical simulations becoming 
available in recent decades. A much richer physics of plas- 
ticity can, however, be expected in magnetoelastic amor- 
phous materials like metallic glasses, and the aim of this 
letter is to commence a research program in this direction. 

The fundamental physics of plastic instabilities in non- 
magnetic amorphous solids has been uncovered in recent 
years in the context of athermal, quasistatic mechanical 
(AQS) strain. While many experiments are done at finite 
temperature and finite strain rates, AQS studies afford 
a unique laboratory for exposing the clean fundamental 
physics of plastic instabilities. In particular bifurcation 
theory, sometime known as catastrophe theory, provides a 
powerful framework for the discussion of the singular be- 
havior of such instabilities. For purely mechanical strains 
irreversible plastic events occur in AQS conditions when 



an eigenvalue of the hessian matrix H. hits zero. The Hes- 
sian matrix is defined as 



n, 



dU 
dridr^ 



Nd x Nd symmetric matrix (1) 



where U({ri(-y)}, 7) is the total potential energy of a solid 
containing N particles whose coordinates depend on the 
external strain 7, i.e. {1*1(7) ■ • • r Jv(7)}- The generic 
plastic event JTJ[3] is a codimension-1 saddle-node bifur- 
cation which occurs at some value of 7 = 7p of the ex- 
ternal strain where the lowest eigenvalue of the Hessian 
matrix (excluding Goldstone modes when they exist) hits 
zero like s/jp — 7. The saddle-node bifurcation has been 
shown to be universal for a large variety of amorphous 
solids from simple binary glasses with pair potentials to 
metallic glasses with many-body interaction potentials [3] . 
Among other interesting phenomena this square-root sin- 
gularity determines the system-size dependence of the av- 
erage stress drops (Act) and average energy drops (AC/) 
during the plastic events in the elasto-plastic steady state 
of amorphous solids above the yield stress [1] , 



(Act) ~ N p , (At/) ~ N a 



(2) 



with universal exponents a = 1/3, j3 = —2/3. 

Modeling magneto-plasticity: The aim of this Let- 
ter is to discuss the rich and interesting physics of plastic 
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instability and failure in general in magnetic amorphous 
solids in which two control parameters exist, i.e. the exter- 
nal strain 7 and a magnetic field B. We cannot cover here 
all the exciting physics that involve cross-effects between 
strain 7 and magnetic field B in phenomena such as mag- 
netostriction, the Villari effect, Barkhauscn noise, mag- 
netic domain evolution and consequently their effects on 
the mechanical properties of materials like metallic glasses 

Our aim in this Letter is to point out to the interested 
community that magnetoelastic materials have a lot to 
offer in terms of new riddles, and here we focus on the B— 7 
phase diagrams including the existence of codimension-2 
plastic instabilities. Possible totally new singularities will 
be discussed at the end. 

In amorphous solids whose constituents are magnetic 
the potential energy U consists of two parts 

U({r i ( 1 ,B)},{S i ( 1 ,B)}, 1 ,B) = U mech ({r i ( 1 ,B)}, 1 ) 
+£W({n(7, B)}, {5i( 7 , B)}, 7, B), (3) 

where Si are spin variables and the mechanical part can 
be written as 

EUch({ri(7,B)},7) = E «(ry(7,S),7) . (4) 

where r,j = \ri— Tj\. (in this Letter we shall treat both the 
strain 7 and magnetic field B as scalars) . The interparticle 
potential typically exhibits a minimum at a scale ry = a. 
Note that most glassy materials are made of a number of 
constituents and therefore this scale will be different for 
each distinct pair of constituents. 

The magnetic part needs to be modeled to fit best a par- 
ticular material, and different materials will have some- 
what different magnetic interactions. For concreteness 
therefore let us consider here an amorphous 2-dimcnsional 
model of classical xy spins in which the orientation of each 
spin is given by an angle fa. We also denote by Oifa) the 
local preferred easy axis of anisotropy, and end up with the 
magnetic contribution to the potential energy in the form 
(TTj : where J(rij ) is the exchange interaction overlap in- 
tegral whose rij dependence needs to be found for each 
material. As shown in appendix A, the exchange integral 
is expected to have a relatively sharp peak at a typical 
scale r* . Of crucial importance is the relative dimension- 
less ratio r* J a where a is the position of the minimum 
set by interparticle potential (in a typical glass this ratio 
will be different for every given pair of constituents). Con- 
sider for example a small uniaxial compressive stress on a 
given material. For r*/a > 1 this will result in decreas- 
ing the magnetic interactions with increase of compres- 
sive stress. On the other hand for r* /a < 1 compressive 
stress will result in increasing the magnetic interactions. 
Indeed, in this case it is even possible that larger compres- 
sive stress will turn around this behavior. In the case of 
tensile stresses the oppositive behavior can be expected; 



while for shear stresses anisotropic phenomena might ap- 
pear. 

The second term of Eq. (0 that models the effect of local 
anisotropy will tend to form magnetic domains, depending 
on the relative magnitude of the parameter K with respect 
to Jq where q is the average number of nearest-neighbors 
(note that K could in principle also be a function of r^j , 
but we expect that the existence of local anisotropy due 
to the local positions of the atomic neighborhood is more 
important.) For K/ Jq > 1 we expect a complex inter- 
play between domain formation, a vortex glass phase and 
ferromagnetic behavior. 

The last term is the interaction with the external field 
B. here fa = fasgiLi where fis is the Bohr magneton, 
Li the angular momentum of the state considered and gi 
is the associated Lande g-factor. Note that the first two 
terms will exhibit important coupling between strain and 
magnetism due to the dependence on r^. Note also that 
the macroscopic magnetization M of our sample is M = 

v ^->i l J " i ' 

At this point the joint dependence of the potential U 
on the positions rj and the angles fa clarifies; it implies 
that strain 7 and magnetic field B can be changed quasi- 
statically to maintain the magnetic amorphous solid in 
mechanical equilibrium. This requires that the net forces 
and torques on each particle vanish , 



dU_ _ dU_ 

dr., dfa 



(6) 



for all strains and magnetic fields. Thus, if the system 
is subject to an affine transformation r, — > ^(7) ■ r,, the 
particles will exhibit additional non-affinc displacements 
to annul the resulting forces, but then also the spin ori- 
entations will change to minimize the energy, and in total 
we expect a non affine response according to 



h(i) ■ r % + Ui^ NA 



,NA 



(7) 



Similarly, in response to a change in the applied magnetic 
field B both the particles positions and the spin orien- 
tation will change; we denote the total changes that are 
required to keep the system in mechanical equilibrium as 
u hNB and faNB' 



r t + Ui, NB 



+ < 



H.NB ■ 



(8) 



Plastic instabilities: To proceed further, we need to 
study the generalized Hessian matrix "H, which is now of 

rankiV(d+l) x N(d+l): 



n 



n (rr) 
n (4>r) 



■Tj(<t"t>) 



(9) 
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U mag ({rj},{Si}) = - y^ J(r\j) cos I 



"» Vj 



- K y^ cos (<fo - 6>i({rj})) 2 --By^jUjCos^) . 



<u> 



(5) 



Here the four partial matrices are defined by 

a 2 u 



n 



(rr) 



u 



{r<t>) 



u 



(4>r) _ 



n 



(4>4>) 



dridrj 
d 2 U 


dr l d(j) J 

d 2 u 


d&drj 
d 2 U 



(dN x dN matrix) 
(dN x N matrix) 
(N x dN matrix) 
(N x N matrix) . 



(10) 



The Hessian matrix 7i is real and symmetric, and there- 
fore diagonalizable. Besides Goldstone modes with zero 
eigenvalues all the other modes are associated with posi- 
tive eigenvalues as long as the system is mechanically sta- 
ble. Plastic instabilities will occur when the lowest posi- 
tive eigenvalue of W will approach zero upon the increase 
of 7, B, or both. These instabilities are manifest in the 
equations of motion of the non-afhnc responses defined in 
cqs.© and ©. These are denoted as 



,(7.' 



.0>,r) 



# U J, 



NA , 



,(7,0) 



do 



'3,NA | 



c*7 



Qj 

(6,0) _ d0j,NB | 



(11) 



j dB n ' o as ' ' 

Note the physical meaning of these non-afhne flows. The 
first is the well known non-affine response which occurs 
also in non-magnetic amorphous solids as a result of an 
external strain. There is now a new response to the ex- 
ternal strain, namely that of the spins. Finally we have 
two new non-afflnc responses to the magnetic field, i.e the 
response of the spins and the response of the particle po- 
sitions. One can derive the equation of motion for these 
flows in the form 

v (7) = fi 1 ■ S^ 



Ah) 



1 , w (b) 



-n-'s 



(12) 



where the mismatch forces are defined in Appendix B and 



,(7) 



■1,(7,0) 



Ah) 



(b,r) 



(13) 



It is now manifest in Eq. (|12[) that when an eigenvalue 
of H hits zero the non-affine response can be very large 
indeed. 

The derivation of the equations of motion of the eigen- 
values of the Hessian matrix follows verbatim the tech- 
niques uses in the purely mechanical case |12j . These read 

^fftft + ^l 



<9 7 ' 

dX k 

dB 1 



.(7) 
-kk 



kk 



H a " 



J kkl 



'kkf'i 



E 



M 



[b 



(r) 

kkf. ' 



u kkl\ 



(14) 



The coefficients in these equations are all defined in Ap- 
pendix B. 

Symmetries and Codimension-2 Instabilities: In 

the rest of this Letter we focus on the novel co-dimension 
2 plastic instabilities that arise in the present model. To 
this aim note that while the potential energy U has no 
special behavior under 7—^—7 for each configuration of 
{r;},{0i}, it does remain invariant under the combined 
transformations 



[{n}, {&}; 7, B] -> [{n}, {& + tt}; 7, -B]. 



(15) 



This symmetry has important implication for the proper- 
ties of strained amorphous solids in the presence of mag- 
netic fields. One immediate consequence is that the eigen- 
values of the Hessian matrix obey this symmetry as well, 

A fe ({r i },-M;7,--B) = A fe ({r- i },M;7,B), (16) 

where we have now returned the vector notation to the 
magnetic field. Consider now the B,j parameter space, 
firstly along the line B = 0. Along this line we expect 
to find a first plastic instability at a value of 7 = jp. 
It is easy to see from Eq. [H] that for B = and when 
one of the eigenvalues tends to zero, say Ap — » 0, the 
solution has the square-root singularity, Xp ~ VTp — 7, 
which is the hallmark of a saddle node bifurcation. From 
this point there emanates a line of instabilities at values 
of 7 = r )p{B) where (again from examining Eq. (|14l0 . 



Ap(7,B)^ v / 7 p(B)-7. (17) 

Due to the symmetry (|16[) we expect this line to depend 
on the value of the magnetic field according to 

7p(B) = 7p - fciM ■ B - k 2 B 2 + ... ifM^O 
lP (B) = 7p - k 2 B 2 - k A B A + ... if M = 0. (18) 

A consequence of Eq. (fTg)) is that for a small magnetic 
field B the lowest plastic mode will become unstable at a 
finite magnetic field 

Bp{l) = (7p-7)/(*iA0 ifM^O 
Bp{i) = y/(j P - i)/k 2 ifM = 0. (19) 

for 7 < 7 P but close to it. Inserting Eq. (fH?j) in Eq. ([T5| 
and then in (|17[). and then expanding to first order, we 
find 



Ap( 7 ,B) - Vm^(B p (j) - B) ifM^O (20) 

\p{l,B) ~ (7p-7) 1/ V(^p(7)-- b ) ifM = 0. 

At this point we examine a co-dimension 2 path in the 7, B 
plane, denoted as B(j). Along this path the eigenvalue 
satisfies the equation 

dXp/dj = (-^-\ B ) + (-^)| 7 (dB( 7 )/d 7 ). (21) 
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Now if the curve .6(7) is analytic the saddle- node expo- 
nents will not change, but consider the singular curve in 
the (7, B) plane that ends at (jp, Bp) of the form 



(B P -B) = K(jp(B P ) - j) 6 



(22) 



Now (dB/dj) = ~5K{ lP {B P ) - ^f- 1 and different be- 
havior can be expected for S < 1 and 8 >= 1. If 6 >= 1 
the saddle-node exponents remain unchanged but if 6 < 1 
then 

Ap~( 7P (Bp)-7) 5 / 2 . (23) 

This is the main result of the present calculation, showing 
that co-dimension 2 paths can yield a range of singularities 
to the lowest eigenvalue provided S < 1 . As a consequence 
the exponent characterizing the approach of Xp to zero 
can now in principle take any value smaller than a 1/2. 

The discussion so far has been based on two implicit as- 
sumptions that we want to discuss here further. The first 
assumption concerns the thermodynamic limit TV — > 00. 
Because we are dealing with amorphous solids we might 
expect that even in the case that M = in the thermo- 
dynamic limit, for finite sample M ~ 1/V7V. Thus in the 
presence of small magnetic fields the first of Eqs. (18) 
is always dominant. The other assumption was that the 
first plastic instability at 7p is well separated from a sec- 
ond plastic instability at, say, 7p. When B increases the 
instability lines emanating from these two points at B = 
can come closer and maybe even cross. This is definitely 
an interesting bifurcation where two eigenvalues go to zero 
simultaneously, bringing about the need for a new theory 
and new consideration [13j . 

The existence of a range of scaling behavior for the 
eigenvalues of the Hessian that hit zero has important 
implications on the statistical physics of magnetoplastic 
materials, putting them in a separate universality class 
from the strictly mechanical instabilities seen in elasto- 
plastic materials. We note that we only discussed the case 
of weak magnetic fields where the nonlinear terms in Eq. 
(fTHl) are negligible. Note that in these equations the sign of 
the various coefficients fci , ki and k± is crucially important 
in determining different possible phenomcnologics. Second 
order perturbation theory suggests that k 2 > |T3]. De- 
pending on the signs of these coefficient there may exist a 
point (pf c ,B c ) at which this line terminates |13j . In such 
a case we expect to see very interesting scaling behavior 
whose precise nature must await future work. 

Appendix A: The typical exchange integral. 

For solids where only small deviations from the equilib- 
rium positions occur this integral can be given by the two 
first terms of a Taylor expansion around the equilibrium 
atomic spacing tq at zero temperature (see, e.g., P3HHII) 



J(r) = J + J'(r-r ), 



(24) 



where the constant J' is responsible for coupling between 
clastic and magnetic energies and can be positive or neg- 
ative. For amorphous solids in which the particles may 



migrate relatively large distances this is not a sufficient 
model. The simplest model that appears consistent with 
observation has an exponential long-range decay (see, e.g., 
|17U18j ) and a rapid rise at short distance, for example 



J(r) = J e- /3(r - ri) , r>n 
J{r) = J e Q(r - ro) , r<r . 



(25) 



The rise at short distances together with the decay at 
large distances result in a peak at intermediate distances 
in agreement with observations |19j . In the intermediate 
range the exchange integral can be fitted by a polynomial 
to ensure continuity and differentiability. 

An example of the exchange integral function for Nickel 
based on the Monte Carlo / Molecular Dynamics results of 
[IB] is shown in Fig. Q] Note that J{r) is not a monoton- 
ically decreasing function of r. In fact, the overlap inte- 
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Fig. 1: The exchange integral computed by fitting to experi- 
mental data in Nickel. The interparticle distance is measured 
in units of a which is the length appearing in the mechanical 
two-body interaction. The parameters a and /9 were defined in 
Eq. (A2). 

gral tends to peak roughly at the typical distance between 
particles, set by the mechanical interaction (fj|, denoted 
as a in the ordinate of Fig. Q] Thus the effects of strain 
and magnetic fields will be crucially dependent on where 
exactly J(r) peaks with respect to a. For example the 
present form of J(r) means that for crystalline Nickel the 
magnetization should decrease under increasing pressure. 

Appendix B: Calculation of Coefficients. The 

coefficients in Eqs. Q3]are determined by 

where ?/r ' is the eigenfunction of the Hessian associated 
with eigenvalue Afe, 



;(7) = 


/w(7,>-)\ 
( S (7>0) ) ' 


- ~ I S (W) J 


(7,0 


d 2 u 

djdr 7 


-(7,0) - d2{J 


<b,r) 


d 2 U 


o(b,<j>) _ d U 



(27) 



dBdr 



dBdcj) 



(28) 
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a (7) = B ( 7 ) . ^k) c (7) = ^ : ^M ; W = S (6) . ^ W ^ C W = H : ^W 1 6 W n = TKr) : ^W^m) 

bil = r<<«:v> { W m) , (26) 



Finally, [19] N. Sakai and H. Pifer, Phys. Rev. B 33, 1875 (1986). 



I \<r{<t>rr) <r{<t><t>r) J ' ' \q-(4>r<t>) j-(4>4>4>) J 

(29) 
and we have denoted the third-order tensor of potential 
energy derivatives as, for example, 7 ~S = dr g^.g r an d 

' ijk — dfadfadn- 
REFERENCES 



[1] D. L. Malandro and D. J. Lacks, Phys. Rev. Lett. 81, 

5576 (1998). 
[2] C. Maloney and A. Lematre, Phys. Rev. Lett. 93, 016001; 

93, 195501 (2004). 
[3] R. Dasgupta, S. Karmakar and I. Procaccia, Phys. Rev. 

Lett. 108, 075701 (2012). 
[4] S. Karmakar, E. Lerner, and I. Procaccia, Phys Rev. E 

82, 055103R (2010). 
[5] R.L. Steinberger, Physics 4, 153 (1933). 
[6] J. Kamarad, Z. Arnold, J. Schneider and S. Krpicka, J. 
Magnetism and Magnetic Materials 15-18, 1409 (1980). 
[7] J.D. Livingston, Phys. Stat. Sol. (a) 70, 591 (1982). 
[8] D.C. Giles, J.Phys. D, Appl. Phys. 28, 1537 (1995). 
[9] G. During and S. Zapperi, J. Appl. Phys. 85, 5195 (1999) 
[10] L. Lanotte, G. Ausanio, M. Carbuchiccio,V. Iannotti and 
M. Muller, J. Magnetism and Magnetic Materials 15-18, 
276 (2000). 
[11] The savvy reader might wonder where is the long 
range dipole-dipole interaction energy Udd({i'i},{4'i}) = 
-5Z <ii> [3m i • lYjirij • ry - m; • mjTyJ/rjj,-, which will 
encourage magnetic domain formation in the amorphous 
solid. It is clearly of vital importance for domain forma- 
tion in crystalline solids and in materials made of grains 
of super-paramagnetic particles, but in the case of amor- 
phous solids is likely to be less important than the local 
anisotropy term and will be neglected here. 
[12] H.G.E. Hentschel, S. Karmakar, E. Lerner and I. Procac- 
cia, Phys. Rev. E 83, 061101 (2011). 
[13] H.G.E Hentschel and I. Procaccia, "Theory of magneto- 
plastic materials", in preparation. 
[14] Z. Y. Chen, M. Kardar, Elastic antiferromagnetics on a 

triangular lattice. J. Phys. C 19, 6825-6831 (1986). 
[15] L. Gu, B. Chakraborty, P. L. Garrido, M. Phani, J. L. 

Lebowitz, Phys. Rev. B 53, 11985-11992 (1996). 
[16] B. Grossmann, D. G. Rancourt, Phys. Rev. B 54, 12294- 

12301 (1996). 
[17] D. K. Belashchenko, V. V. Hoang, P. K. Hoang, J. Non- 

Crystall. Solids 276, 169-180 (2000). 
[18] R. E. Coffman, G. R. Buettner, . J. Phys. Chem. 83, 
2387-2392 (1979). 



